The Astrophysical Journal, in press 

Preprint typeset using L^T^X style cmulatcapj v. 5/14/03 



RADIATIVELY INEFFICIENT MHD ACCRETION-EJECTION STRUCTURES 
Fabien Casse and Rony Keppens 

FOM-Institute for Plasma Physics Rijnhuizen, Association Euratom/FOM, 
P.O. Box 1207, 3430 BE Nieuwegein, The Netherlands 
fcasse@rijnh.nl, keppens@rijnh.nl 
Submitted, July 18th, 2003; Accepted October 2nd, 2003 
The Astrophysical Journal, in press 

ABSTRACT 

We present magnetohydrodynamic simulations of a resistive accretion disk continuously launching 
transmagnetosonic, collimated jets. We time-evolve the full set of magnetohydrodynamic equations, 
but neglect radiative losses in the energetics (radiatively inefficient). Our calculations demonstrate 
that a jet is self-consistently produced by the interaction of an accretion disk with an open, initially 
bent large-scale magnetic field. A constant fraction of heated disk material is launched in the inner 
equipartition disk regions, leading to the formation of a hot corona and a bright collimated, super- 
fastmagnetosonic jet. We illustrate the complete dynamics of the "hot" near steady-state outflow 
(where thermal pressure ~ magnetic pressure) by showing force balance, energy budget and current 
circuits. The evolution to this near stationary state is analyzed in terms of the temporal variation of 
energy fluxes controlling the energetics of the accretion disk. We find that unlike advection-dominated 
accretion flow, the energy released by accretion is mainly sent into the jet rather than transformed into 
disk enthalpy. These magnetized, radiatively inefficient accretion-ejection structures can account for 
under- luminous thin disks supporting bright fast collimated jets as seen in many systems displaying 
jets (for instance M87). 

Subject headings: Accretion, accretion disks — galaxies: jets — ISM:jets and outflows — MHD 



1. ACCRETION DISKS AND JETS 
1.1. Accretion- Ejection models 

Astrophysical jets are quite common phenomena across 
our visible universe. They are typically observed in as- 
sociation with accreting objects such as low-mass young 
stellar objects (YSO), X-r ay bin a ries ( XRB) or active 
galactic nuclei (AGN, e.g. iLiviol l|1997l) and references 
therein). In all these systems, the mass outflows exhibit 
very good collimation at large distances from the central 
object as well as high velocities along the jet axis. Al- 
though operative at widely disparate lenghtscalcs, these 
accretion-ejection systems have other features in com- 
mon, in particular, observational links have been es- 
tablished in all cases between accretio n disk luminos- 
ity an d jet emission: for YSO see e.g. [Hartigan et aP 
IU991. for XRB see lMirabel et alJ (|1998j) and for AGN 
ISerieant et aP ljl998j) . 

The most promising unifying model relies on a sce- 
nario where an accretion disk interacts with a large- 
scale magnetic field in order to give birth to bipolar 
collimate d jets (for t he spe cific case of early protostars 
se e also iLerv et alJ (119991)^1. Since the seminal work 
bv lBlandforclfcF'avnel l)1982|) . it is known that the action 
of an open magnetic field configuration threading a disk 
can brake rotating matter in order to transfer angular 
momentum into the jet and provide energy for accelera- 
tion of jet matter. This magnetohydrodynamic (MHD) 
model describes the interaction of the accretion flow with 
a magnetic field whose ori gin can be due to a dvection of 
interstellar magn etic field llMouschoviasll976t) and/or dy- 
namo produced l)R.ekowski et alJl2003l) . The collimation 
of the flow is self-consistently achie ved by the electric cur- 
rent produced by the flow itself l)Hevvaerts fc Normanl 
1989). This current provokes a radial pinching of the 



plasma that can b alance both magn etic and thermal 
pressure gradients iSautv et alJ I 2002D . The analytical 
work bv lBlandford fc Pavnei l)198 2) did make the simpli- 
fying assumption of a cold plasma. 

Numerous studies have dealt with these Magnetized 
Accretion-Ejection Structures (MAES) in the last two 
decades. Sophisticated semi-analytical models deal with 
stat ionary self-similar investiga tions gradually extending 
the Blandf ord fc Pavnei l|1982l) model with more physi- 
cal effects, starting from simple v ertical mass flux pre- 
scriptions (|Wardle fc Konig]|ll993ft to realist ic disk equi- 
libria where a mbipolar diffusion (0 11996(1 , resistivity 
l)Ferreiralll997fl or both viscosity and resistivit y are ad- 
equately incorporated l|Ca,sse fc Ferreiral EKKl . These 
studies bring deep insight in the physical conditions 
prevailing at the disk surface required to launch jets, 
but fail to give realistic jet topologies in the super- 
Alfvenic region. The latter has spurred a variety of 
numerical MHD studies, aiming at a more realistic de- 
scription of the trans- AlfVenic flows. The complex- 
ity of the dynamics of the accretion-ejection flow has 
forced many authors to either focus on jet dynamics 
alone llKrasnonolski et al.lll999l lOuved fc PudritzllT997l 
lUstvugova et alJ 119951 11999ft or to s tudy disk-outflow 
dynamics over very short timescales llKa to et al l l2002t 
iMatsumoto et alJ 119961: lUshida fc Shibatal H985fT The 
latter class of studies were done using ideal MHD frame- 
work which is inconsistent with long-term jet produc- 
tion since the frozen-in magnetic structure is advected 
with the disk material. This leads to a rapid destabi- 
lization of the system due to a magnetic flux accumula- 
tion in the inner part of the disk. Moreover, in most of 
the previous studies, the accretion disk is modeled as a 
non-Keplerian thick torus without any outer mass inflow 
which would mimic the mass reservoir of the disk outer 
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regions. Here again this compromises the long-term jet 
production since accretion will end rapidly because of the 
lack of disk material. 

In contrast to these studies, iCasse fc Keppensl l|2002|) 
(hereafter CK02) recently presented MHD computations 
demonstrating the launching of super- fastmagnetosonic, 
collimated jets from a resistive accretion disk over a large 
number of dynamical timescales. These axisymmetric 
simulations heavily relied on the general requirements 
identified by self-similar analytical models, in particular 
the necessity for equipartition disk r egions together with 
suffic i ently bent magnetic surfa ces fBlandfor d & Pavnel 
119821 iFerreira fc Pelletieilll995|h Nevertheless, for sim- 
plicity CK02 replaced the energy equation by an adi- 
abatic polytropic relation which did not enable us to 
explore the energetics of the resulting accretion-ejection 
flow. This is of primary interest for realistic modeling, 
since accretion disks supporting jets display a disk lumi- 
nosity much smaller than expected in the framework of 
any currently accepted 'standard' model. 

1.2. Under-luminous accretion disks and hot collimated 

jets 

The detection of accretion-ejection motions is typi- 
cally achieved through multi- wavelength observations. 
Emissions from the disk and from an associated out- 
flow differ since they are produced by different mecha- 
nisms in media of widely different densities and tempera- 
tures. Accretion disks displaying outflows always present 
low-radiative efficie ncy, and this is true in binary stars 
(IRutten et al.lll 992h microquasars l|Mirabel et alJH998j) 
or AGN IpiMatteo et all 120031 12000|) . One scenario to 
explain the lack of disk emission invokes a weak cou- 
pling between ions and electrons, resulting in electron 
temperatures being much lower than ion temperatures 
()Rees et alJ Il982t iShapiro et alJ fl97^ . This scenario 
has led to the concept of "a dvection-dominated a ccre- 
tion flows" (ADAF, see e.g. iNaravan fcTil l| 19951) and 
references therein) where energy released by accretion of 
matter is primarily transferred to heat the ions which 
are ultimately advected onto the central object. If this 
is a black hole, this energy is t hen simply lost. Recently, 
iBlandford fc Beeelmanl l)1999|) presented an alternative 
to ADAF, where inclusion of an outflow from the disk 
was suggested, in which case a part of the accretion en- 
ergy is sent into outflowing hot matter. Although several 
numerical MHD simulations of radiatively inefficient ac- 
cret ion disks have now been reported in the liter ature 
(see lleumenshchev. Naravan k, Abramowicd l)2003|) and 
references therein), they have so far been unable to de- 
scribe self-consistently both radiatively inefficient accre- 
tion flow and persistent, hot and collimated jet outflow. 

The aim of the present paper is to compute radiatively 
inefficient magnetized accretion-ejection flow where both 
the inward accretion and the upward, transmagnetosonic 
ejection flow can provide a way out for the energy re- 
leased by accretion. This is achieved by fully accounting 
for an appropriate energy equation in the MHD descrip- 
tion of the interaction of an accretion disk with a large 
scale magnetic field. We present the theoretical and nu- 
merical background of our model in Section 2. In Sec- 
tion 3, we display results obtained from simulations of 
a resistive accretion disk with full energy consideration. 
In Section 4, we analyze the temporal evolution of the 



energy balance including both disk and jet powers. We 
conclude in Section 5 with a summary and open issues 
to address in future work. 

2. MAGNETIZED ACCRETION-EJECTION FLOW 

We aim to model a magnetized accretion disk launch- 
ing jets. The description of its plasma dynamics is 
achieved in an MHD framework. We first present the 
full set of resistive MHD equations, then discuss initial 
and boundary conditions. We conclude this section by 
explaining the physical normalization employed in our 
simulations, since this makes them applicable to very dif- 
ferent systems (YSO to AGN). 

2.1. MHD equations 

The full set of MHD equations expresses conservation 
of macroscopic quantities as mass, momentum and en- 
ergy. The mass conservation is 
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where p is the plasma density and v the velocity of the 
fluid. Since we are assuming axial symmetry, note that 
only the poloi'dal component of velocity is involved in 
Eq. Q . The momentum conservation accounts for forces 
arising from thermal pressure P, Lorentz force, centrifu- 
gal force and gravity, namely 

^ + V-(vpv-BB) + v(^ + p)+pV$ G = 0, (2) 

where B is the magnetic field and $g = —GM*/(R 2 + 
Z 2 ) 1 / 2 is the Newtonian gravity potential. Coordinates 
(R, Z) stand for cartesian coordinates in the poloi'dal 
plane. The temporal evolution of the magnetic field is 
governed by the induction equation which must take into 
account the effect of a temporally and spatially varying 
resistivity 77, 
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where J is the current density defined by the Maxwell- 
Ampere equation 

J = V x B . (4) 

We have adopted units where vacuum permeability p is 
equal to unity. The total energy density e is defined as 
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where 7 = 5/3 is the specific heats ratio of a non- 
relativistic plasma. The energy equation in a resistive 
framework then reads 



de 
W 



P 



2 



BB-v =? 7 J 2 -B-(VX7 ? J). 

Note that resistivity, operational in the disk alone, pro- 
vides a local Joule heating. Ohmically heated plasma can 
then be advected by both accreting and ejecting flow in 
the case of a radiatively inefficient plasma. The last rela- 
tion closing the set of equations links pressure and den- 
sity to temperature T = Pj p, as we assume the plasma 
to be ideal. 
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2.2. Initial conditions 

The initial configuration of the system is very close to 
the one adopted in CK02. The computational domain 
is [R, Z] = [0, 40] x [0, 80] with a resolution of 104 x 204 
cells and includes a sink region around the origin to avoid 
the gravitational singularity. The size of the sink region 
in CK02 was several grid cells long (up to the inner disk 
radius setting our unit length) and only two cells high. If 
one wants to simulate systems ranging from YSO (New- 
tonian gravity) to b lack holes (with for instance a pseudo- 
Newtonian gravity IPaczvnskv fc Wiital l)198fl) . suitably 
describing plasma dynamics from several Schwarzschild 
radii outwards), we have to take into account the fact 
that the gravitational singularity does not always occur 
exactly at the origin but can be located at some distance 
from it (in the case of a black hole). Therefore, we now 
modify the size of the sink region, and design it to cover 
one unit length both along R and Z . For a complete de- 
scription of the boundary conditions, for both the treat- 
ment of this interior sink region as well as the domain 
boundaries, we refer to CK02. 

We briefly list our nearly unchanged initial conditions 
for density and velocity, namely 
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Note that in the present paper we have decreased the 
initial density range compared to CK02 for computa- 
tional performance. However, we still cover 6 decades 
in density contrast at t = 0. Similar to CK02, we as- 
sume the initial thermal pressure to be polytropic such 
that P = p 5 / 3 . The disk height H = eR increases lin- 
early with radius, with aspect ratio e = 0.1 suitable for 
a thin disk. The constant R Q — 4, while the parame- 
ter m s = 0.1 is consistent with initial subsonic accre- 
tion motion. The ms parameter range is restricted by 
both subsonic accretion motion and the minimal bend- 
ing of the polo'i dal magnetic surface in or der to achieve 
jet production l|Blandford fc Pavnel 11982). Indeed the 
induction equation shows that a steady-state magnetic 
configuration with bent magnetic surface above a thin 
disk requires both resistivity and minimal radial motion 
since 

dB 

V-gf ~ -VrB z . (8) 

The choice of m s is then to be smaller but close to unity. 
The initial magnetic field configuration is purposely 
taken very different from CK02. Where CK02 started 
from a radially stratified, but purely vertical field, here 
we start with an open polo'idal magnetic structure to ad- 
dress the issue of jet collimation adequately. The main 
difficulty is to find an initially bent magnetic configura- 
tion which reconciles several constraints like symmetry 
at the equatorial disk plane and at the jet axis with the 



obvious requirement of V ■ B = 0. Moreover, if one wants 
to have an almost constant plasma beta j3p = B 2 /P ~ 1 
inside the accretion disk, a necessary condition for jet 
launching, the radial variation of the magnetic compo- 
nents must be close to a power-law of index —5/4. In 
order to fulfill the above statements we adopt the mag- 
netic configuration ((3p = 0.6 in our simulation) which 
exploits the function 



F(R, Z) . 



and given in full by 
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The extra term appearing in the Bz definition avoids 
vanishing vertical magnetic field at the jet axis (this term 
rapidly vanishes as R increases) and does not interfere 
with the constraint V • B = 0. The initial magnetic 
field configuration is not force-free but prevents from any 
initially-induced collimation since magnetic surfaces tend 
to widen (outwardly decreasing magnetic pressure). On 
figure |QJ, we have displayed this initial polo'idal mag- 
netic configuration, where we have set C = 0-04 (£ con- 
trolling the initial bending of magnetic surfaces). 

The total energy density e is then found from the previ- 
ous quantities according to equation JSJ. The anomalous 
resistivity r\ is prescribed in the same way than in CK02, 
namely 



n = a m V A \ Z =Q Hexp 



H 2 J 
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which is an a-prescription (a m = 0. 1) where the Alfven 
veloc ity replaces the sonic speed ijShakura fc Sunvaevl 
|1973|K Through the dependence on the Alfven velocity, 
this becomes a spatio-temporally varying profile which 
essentially vanishes outside the disk. 

2.3. Numerical code: VAC 

The numerical calculations presented here are done us- 
ing the Versatile Advection Code (VAC. see Toth (1996) 
and http : //www.phys .uu.nl/ ~toth). We solve the full 
set of resistive MHD equations under the assumption 
of a cylindrical symmetry. The initial conditions de- 
scribed above are time advanced using the conservative, 
second ord er accurate Total Varia tion Diminishing Lax- 
Friedrichs ijToth fc Ods trcil 1996) scheme with minmod 
limiting applied on the primitive variables. We use a 
dimensionally unsplit, explicit predictor-corrector time 
marching. To enforce the solenoidal character of the 
magnetic field, we apply a projection sc heme prior to 
every time step l)Brackbill fc Barneslll98fl|) . 

2.4. Accretion disk properties 

The initial conditions discussed in section 12.21 are in- 
spired from self-similar descriptions of accretion disks. 
It is important to realize that through the careful choice 
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of normalization used in these expressions, our results 
will be applicable to both YSO and AGN systems. In- 
deed, units are now expressed in terms of inner disk radii 
R = 1, velocities from QkH |h=i, and density values 
(very close to or) at R = 1. We here list what this di- 
mensionalization implies for both black hole and YSO 
systems. 

In this particular f ramework, the standa r d ac- 
cretion disk mod el lILynden-Bell fc Pringle) 11974 
iNovikov fc Thornd 11973 iShakura fc Sunvaevl Il973|) 

describes astrophysical accretion disks provided that 
the accretion is sub-Eddington. If radiative pressure is 
neglected, it can be shown that the plasma temperature 
scales as 
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where flm p is the mean molecular mass of the plasma 
(proton mass m p ), H = eR is the disk height, il^ = 
(GM^/i? 3 ) 1 / 2 the Keple rian angular velocity and ks 
the Boltzmann's constant. The first expression contains 
Rs = 2GM*/c 2 , which stands for the Schwarzschild ra- 
dius, and is applicable in the case of a black hole. The 
second expression is to be used for a YSO disk. 

This temperature scaling combined with the definition 
of the accretion rate M a leads to all other dimensional 
quantities. When applying our results to a certain sys- 
tem with inner radius R, central object mass M» and 
accretion rate M a , the dimensional scale factors are: 
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Again, the first expressions are to be used for black hole 
systems, while the second lines are for YSO. 



The last quantity is the magnetic field that we set to 
be close to equipartition with thermal pressure, a neces- 
sary condition for jet launching. Indeed, the above tem- 
perature statement implies that thermal pressure natu- 
rally overcomes the gravitational pinching in standard 
accretion disk (P ~ p£l 2 K H 2 ). Magnetized disks are also 
prone to magnetic pinching so a natural equilibrium of 
the disk is such that P > B 2 / ' \i Q . The former constraint 
coupled with the MHD Poynting flux amplitude thresh- 
old to produce super- fastmagnetosonic, collimated out- 
flows leads to a disk con figuration where P ~ B 2 //.i 
l|Ferreira fc Pellet ierlll995|) . The magnetic field will then 
be 
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Observations indicate that the inner parts of AGN and 
YSO disks mainly differ in their temperature (AGN 
> lO 9 ^, YSO < lO 4 ^). This difference is also evi- 
dent in other properties (pressure, magnetic field), but 
its main consequence is that these systems have very dif- 
ferent dominant radiative mechanisms operational, since 
the opacity regime varies from fully to partially ionized 
plasma. 

3. ACCRETION DISKS AND "HOT" JETS 

In this section, we present the MHD simulations result- 
ing from the initial configuration presented in the previ- 
ous section. We first describe the temporal evolution of 
the accretion-ejection structure, then display the full jet 
launching mechanism and finally check the radial force 
balance that is reached in the jet. 

3.1. Temporal evolution of the flow 

The result of numerically simulating the full MHD dy- 
namics of the resistive magnetized accretion disk is dis- 
played in Fig. ^ The six snapshots represent six poloi'dal 
cross-sections of the structure at different stages of its 
temporal evolution. The initial accretion disk configu- 
ration is close to a radial equilibrium where Keplerian 
rotation balances gravity of the central mass. The rota- 
tion of matter twists the initially purely poloi'dal mag- 
netic field lines such that outflows appear at the surface 
of the disk. The outflow is continuously emitted during 
the further evolution of the system. As seen in the den- 
sity levels of the various snapshots, this outflow is less 
dense 0(1O~ 3 ) than the accretion disk 0(1) but denser 
than the surrounding medium 0(1O -6 ) by several orders 
of magnitude. The evolution of both the mass outflow 
and the magnetic field topology gets closer to an equilib- 
rium as time increases since the poloidal velocity com- 
ponent of the outflow becomes parallel to the poloi'dal 
magnetic field. Moreover, the outflow is radially col- 
limated and crosses the Alfven and fast magnetosonic 
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Fig. 1. — Temporal evolution of a resistive accretion disk threaded by a poloi'dal magnetic field. Color levels represent density level while 
solid lines stand for poloi'dal magnetic field lines. The time unit labeling each snapshot is the number of rotations of the inner radius. After 
a few rotations, outflows are escaping from the disk and remain focused despite an initial bent poloi'dal magnetic configuration. Once jets 
are launched, the structure varies little during the remainder of the whole computation. 



Please download the best quality 
version of this paper at 
http://www-laog.obs.ujf-grenoble.fr/~fcasse/MS56638.pdf 



Fig. 2. — Three-dimensional impression of the final near-stationary endstate reached in our simulations. Color levels represent a volume 
rendering of plasma density and translucent surface stands for a magnetic surface anchored at R = 3 in the disk. Yellow and blue lines 
stand for magnetic field lines and flow streamline respectively. As can be seen, magnetic field lines are twisted by disk rotation which 
provokes mass acceleration as shown by the flow streamline that is initially accreting toward the central object and then turns into the jet. 



speeds, to become superfastmagnetosonic well before it 
reaches the top boundary of the computational domain 
(see also CK02). It can then safely be qualified as a jet. 
In Fig. |21 we present a three-dimensional representation 
of the final stationary endstate reached. We visualized 
a magnetic surface which is anchored at R = 3 which 
nicely shows its initial fanning out as well as its eventual 
collimation. Selected fieldlines are shown in yellow, re- 
vealing their helical character. We also plot a streamline 
(in blue) which shows the path of a plasma parcel that is 



initially accreting within the disk in a spiral fashion, until 
it reaches the jet launch area where it gets propelled into 
the jet. Finally, the figure also contains a volume ren- 
dering of the density in accretion disk-jet system, which 
shows that the density structure at the foot of the jet 
is in fact convex since full collimation is only achieved 
from some vertical distance above the disk. In our simu- 
lations full flow collimation occurs beyond Z — 50, which 
corresponds to typically several tens of disk scale heights. 
In order to quantify the temporal evolution of both ac- 
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cretion and ejection flows, we can have a look at both 
accretion and ejection mass rates. These rates also trace 
the efficiency of the magnetic field to remove angular mo- 
mentum from the disk (to allow accretion) and to acceler- 
ate matter in the jet (which creates a vertical mass flux as 
seen on Fig.^). On Fig. 03 we have displayed these mass 
rates as a function of time, Maj and Ma.e being accre- 
tion rates measured at inner Rj = 1 and external radius 
Re = 40, respectively, while Mjet = 27r J^ E RpVzdR 
is the ejection mass rate measured at a given Z above the 
disk surface. On Fig. the inner accretion rate is nor- 
malized with respect to the value of the outer accretion 
rate which by means of our boundary conditions (fixing 
the mass flux) is set to be constant. The strong increase 
of the inner accretion rate with time is a clear transla- 
tion of the role of the magnetic torque in the disk: as 
the simulation begins, the twist of magnetic field lines 
provoked by disk rotation creates a toroidal component 
for the magnetic field which directly leads to a braking 
magnetic torque. The azimuthal braking of matter re- 
duces the centrifugal force and thus allows matter to be 
attracted by gravity at a higher rate. Nevertheless, as 
seen on Fig. [3] the increase rapidly disappears to let the 
inner accretion rate settle to an almost constant value. 
The existence of this plateau accounts for an accretion 
disk where matter rotation and magnetic tension reach 
an equilibrium state. The ejection rate is slowly rising 
and reaches a final value which is of the order of 22% 
of the inner accretion rate. During the final stages both 
ejection and accretion rates remain almost constant and 
enable us to study in detail the accretion-ejection con- 
nection. The final ratio of inner to outer accretion rates 
is of the order of ten. This value larger than unity arises 
from our choice to impose a constant outer accretion rate 
during the simulation, thus not enabling the structure to 
adjust itself to the outer part of the accretion disks that 
lie outside the computational domain. 

Note that these accretion rates have already settled 
on constant values after about 10 inner orbital periods. 
We have evolved the system further to 30 inner orbital 
periods, and as seen in the snapshots from Fig. ^ the 
overall structure remains rather stationary during that 
time. As the jet is confined to within 20 inner radii, the 
final time T — 30 covers 1.3 rotation period as measured 
at the external foot of the jet, or only 0.13 rotation pe- 
riod of the outer radius at R — 40. However, as will be- 
come clear from the further analysis of the force balances 
achieved in the inner disk regions, this is a truly station- 
ary configuration. Note also that due to the large density 
contrasts, this corresponds to several tens of thousands 
of CFL limited time steps, which make our simulations 
computationally challenging. 

Note that our endstate is in contrast to re- 
cent 3 D simulations of radiativ ely inefficient ac cretio n 
flows ijleumenshchev. Naravan fc Abramowiczl (2003) , 
whose authors claim that we are focusing on a transient 
phase of the flow. However, the quasi-steady endstate ob- 
served in lleumenshchev. Naravan fc Abramowiczl l)2003f) 
is one of dominant magnetic field in the entire computa- 
tional domain. This must be due to the imposed contin- 
uous injection of poloi'dal magnetic flux adopted in that 
study. In contrast, our simulations do not assume an in- 
finite reservoir of magnetic flux, and the magnetic field 
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Fig. 3. — Ejection mass rates and inner accretion rates as a func- 
tion of time. The inner accretion rate, normalized to the constant 
value of the outer accretion rate, strongly increases during the very 
early stages and then remains almost constant during the rest of 
the simulation. The ejection rate increases slowly with time, and 
reaches a final value of the order of 20% of the inner accretion rate. 



settles into a stationary state where inward advection is 
balanced by outward slippage due to the disk resistivity. 
Moreover, in our simulation the level of magnetic pres- 
sure is always close to unity either in the accretion disk 
or in the jet. 

3.2. Force balance and current circuit 

The accretion-ejection zone, as already shown in CK02, 
is the most crucial area when simulating the launching 
of jets. A subtle balance between different forces must 
be obtained in order to achieve a continuous emission of 
matter, as well as the collimation of the resulting out- 
flow. As an example, we display in Fig. ^ the different 
forces in the disk equilibrium at the final time of the sim- 
ulation from Fig.^ In this figure, the poloi'dal forces are 
represented along a given magnetic surface already repre- 
sented in Fig. as well as the magnetic torque, poloi'dal 
velocity and temperature. The vector ep occurring in 
Fig. 2] is defined as Bp/|Bp|. The upper-left plot rep- 
resents the total force applied on the plasma which is 
negative in the disk and becomes positive near the disk 
surface. As already shown in CK02, this configuration 
enables the matter below the disk surface to be pinched 
and to remain in an accretion regime while beyond the 
disk surface, matter is accelerated along magnetic field 
lines, leading to a jet. This balance is obtained from 
competing forces where thermal pressure gradient and 
centrifugal force counteract magnetic and gravitational 
pinching. The crucial point of the accretion-ejection con- 
nection is to accurately capture the sign change of the 
sum of all forces in the numerical stationary state. This 
change of sign can only be achieved if the poloi'dal mag- 
netic force along the magnetic surface changes its sign 
near the disk surface, as shown on Fig. Since this 
force is directly related to the magnetic torque, 



(J x B) • Bp = — (J x B) • B 6 



(15) 



this also implies that the magnetic torque changes its sign 
there, as also shown in Fig.0J This configuration is at the 
core of magneto-centrifugal acceleration since toroidal 
acceleration (or braking) determines poloi'dal accelera- 
tion (or braking) of matter. The poloi'dal velocity field 
shown in Fig. 0] confirms the above statements, since this 
field is consistent with an accretion motion in the disk 
(dominant negative radial velocity) while the vertical ve- 
locity becomes dominant beyond the disk surface. The 
angle between the poloi'dal velocity and poloi'dal mag- 
netic field is displayed on the bottom middle panel and 
again confirms the required accretion-ejection configura- 
tion where flow is perpendicular to the magnetic field in 



Radiatively inefficient MAES 



7 



the disk and becomes essentially parallel to it in the jet 
region. 

The plasma temperature levels are displayed in the left 
panel of Fig. [5] and exhibits an interesting behavior. In 
CK02, we were using a polytropic relation which was 
constraining the temperature to be proportional to p 7_1 . 
This invariably led to a jet where the temperature is very 
low because of the relatively low densities in the jet. As 
we now include the full energy equation (and include 
Ohmic heating), this is no longer happening. 

We can see on Fig.[5]that the disk temperature, which 
was initially smaller than unity at this location, has in- 
creased reflecting the locally operative Joule heating and 
gas dynamics. Rewriting energy equation © by sub- 
tracting momentum and induction equations leads to the 
temperature equation that reads 



1 



dT 

at 



— + v c • VT 



TV- v 



(16) 



7 - 1 V 9t ' J p 

where both Ohmic heating and flow dynamics affect the 
behavior of the temperature. Ohmic heating is occur- 
ring inside the accretion disk and inevitably leads to a 
higher disk temperature than in the initial stage. In the 
accretion disk, the velocity divergence is negative since 
V • v = dV z JpZ + OVr/OR + Vr/R where V R is a slowly 
radially varying negative component. This term will also 
account for a local heating of the plasma, even in the 
case where no Ohmic heating would occur. Assuming 
the structure to be close to a steady-state, one has 



v p ■ VT ~ ( 



P 



TV • v)(7 — 1) > 



(17) 



which is consistent with Fig.[5]where isocontours of tem- 
perature are roughly perpendicular to the streamlines 
in the disk. In the jet region, temperature isocontours 
coincide with streamlines as soon as streamlines have 
turned from an accretion motion (Vr < 0) to an ejection 
motion (Vr > 0). It is noteworthy that at this precise 
location we have strong negative V • v which leads to a 
sudden temperature increase. The resulting "hot" jet is 
a general property of this kind of outflow since it does not 
require local Ohmic heating in the corona. We mention 
that performing the same simulation presented here but 
neglecting Ohmic heating in the energy equation also 
gives birth to jets which have similar temperature than 
the disk, although slightly lower than with Ohmic heating 
included. The relatively small influence of the Ohmic 
heating can be explained in our case by the fact that the 
chosen resistivity is very low since r\ = a m HVA < 10 _2 i?. 

It is possible to clarify another aspect of the accretion- 
ejection mechanism by analyzing the current circuit es- 
tablished in the final structure. Indeed, in an axisym- 
metric framework, it is easy to calculate the current / at 
any location thanks to the Maxwell- Ampere law stating 
that 



J • dS 



B • dl = 2nRB g < 0, 



(18) 



where the contour in the last integral is a circle of radius 
R with its center located at the jet axis. Rewriting the 
poloi'dal magnetic force as 

(JxB).ep = -B 6 ^-e P , (19) 



it is clear that this force is entirely controlled by the lo- 
cal variation of the current along the magnetic field line. 
Looking back at Fig. 0] where the magnetic force pro- 
jected along the magnetic field line is represented, we 
see that the above mentioned change of sign of the force 
has a direct implication for the current, namely that the 
derivative of the current along the magnetic field line 
also changes its sign. Looking now at Fig. |SJ where iso- 
contours of / are displayed, we can see that this leads 
to the formation of closed current loops since the cur- 
rent / start to decrease and then increase. This kind of 
current circuit is representative for a working accretion- 
ejection structure where both magnetic braking in the 
disk and magnetic acceleration in the jet occur. In a 
full meridional cross-section through the disk-jet system, 
this i s reminiscent of a 'butterfly' current circuit ijFerreiral 
1997). 

3.3. Jet collimation 

In CK02, we presented MHD simulations of accretion 
disks launching jets where the resulting outflow was colli- 
mated. In that work, we were able to show that the radial 
force balance in the jet was achieved by magnetic tension 
counteracting magnetic pressure effects (a cold jet). Nev- 
ertheless, it was not a priori obvious that this collimation 
is not affected by our choice of an initial purely verti- 
cal magnetic field configuration. In the present paper, 
we purposely choose a different initial magnetic config- 
uration where the poloi'dal surfaces are strongly bent at 
T = as shown in Fig. ^ Still, the resulting jet is well 
collimated, which is a direct proof that the collimation of 
this flow is arising self-consistently. We display in Fig. d 
the radial balance of the jet for the same snapshot than 
in Fig. 0] at Z = 60. Clearly, the total force is vanish- 
ingly small, proving that there is equilibrium in the jet, 
but the situation is quite different than in CK02. Indeed, 
now the thermal pressure is of the same order than the 
magnetic pressure which leads to a qualitatively different 
radial jet equilibrium. The thermal pressure gradient as 
well as the centrifugal force are now balancing the total 
radial magnetic force. This kind of force balance between 
thermal pressure and magnetic force is another clue to 
the occurrence of a "hot" jet where thermal energy is of 
the same order than the magnetic energy. 
In order to rule out any artificial collimation of the 
jet, we have performed the same simulation with differ- 
ent box size and boundary conditio ns for the magnetic 
field. Since ITlstvugova et alJ (^99), it has been shown 
that axisymmetric magnetic field having open bound- 
aries may induce artificial collimation forces due to im- 
posed vanishing radial current in the ghost cells. The au- 
thors showed that one way to obtain accurate numerical 
MHD flow solutions is to switch from open boundary on 
the toroidal magnetic field to a "force-free" one (where 
B. p ■ V(RBq = 0) and to perform the simulation on a 
square domain in R an d Z. We have followed exactly 
the same approach than lUstvugova et alJ l)1999j) and ob- 
tained a solution very close to the one presented in Fig^ 
For instance, we present in FigElthe density and poloi'dal 
magnetic field structure as well as poloi'dal streamlines of 
the flow at T = 30. The collimation is undoubtedly oc- 
curring since the jet magnetic and velocity field are com- 
pletely vertical and super-fastmagnetosonic as the flow is 
escaping from the domain. The fast-magnetosonic sur- 
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Fig. 4. — Plots of forces projected along a given magnetic surface as a function of the curvilinear abscissa along this magnetic surface 
X(a) at T=30. The footpoint of this surface is at R = 3. From left to right and from top to bottom: total force, thermal pressure 
gradient, gravity, magnetic torque, poloi'dal magnetic force, centrifugal force, radial and vertical velocities, angle between poloi'dal velocity 
and poloi'dal magnetic field and density. The key to achieve an accretion-ejection configuration is to have a magnetic configuration where 
the Lorentz force pinches the disk and accelerates matter in the jet. 



face is shown on Fig. [S] with a dotted-dashed line and 
stands for the location where poloi'dal velocity equals the 
local fast-magnetosonic velocity. 

4. RADIATIVELY INEFFICIENT ACCRETION DISK 
LAUNCHING JETS 

The inclusion of the energy equation in the MHD de- 
scription enables us to investigate the energetics of the 
magnetized accretion-ejection structure. For instance, 
in classical hydrodynamical ADAF, it is shown that the 
energy released by accretion of matter is balanced by 
creation of entropy within the disk. In the present simu- 
lations, the occurrence of a dynamically significant mag- 
netic field opens new ways to transfer and to use accre- 



tion energy. We will begin this section by defining the 
different powers relevant for the disk energetics. Subse- 
quently we will analyze the temporal evolution of these 
quantities in order to identify the predominant way of 
energy transport in the flow. 

4.1. MHD Energy equation 

In Section l2~Tl we presented the MHD equations and 
in particular the energy equation J^J. This equation can 
equivalently be rewritten in a full conservative form as 



R 2 

P + — ) - BB 



v + ?yj x B 



= . 
(20) 
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is defined as 



GO 



tsi 4-0 



20 




20 



Fig. 5. — Isocontours of current I oc RBq. The closed current 
loops are direct evidence of the existence of an accretion-ejection 
structure where magnetic field pinches matter i nside the disk and 
accelerates it beyond the disk surface (see Section [3,2l . This circuit 
is one 'wing' of a full 'butterfly' configuration. 



Using Ohm's law, we can replace the current density J by 
the electromotive field E = B x v + r/J and write the gen- 
eral energy conservation in a form which is independent 
of the MHD regime (either ideal or resistive), 



de 



7 



P 



at 



■V- \F 



MHD 



1 P 

= 



E x B 



= 



(21) 



The last term is the MHD Poynting flux. In an axisym- 
metric framework the two fluxes that control the tempo- 
ral evolution of total energy density are (only the polo'idal 
components of these fluxes are relevant): 

• Accretion energy flux F acc : the dominant com- 
ponent in the disk is the radial one since the 
flow components in the accretion disk obey Vg S> 
I Vr I ^> Wz\- In standard accretion disk models 
( Lvn den-Bell fc Pringld 119741: iNovikov fc Thornel 
11973 IShakura fc Sunvaevll973). entropy is negligi- 
ble mside the disk so that this flux can be reduced 
to F acc ~ pv^$g/2. Note that in the case of a jet, 
the energy flux F acc through the disk surface can 
become significant. 

• MHD Poynting flux: represents the energy carried 
by the magnetic field via angular momentum stored 
in the magnetic field. At the disk surface, where 
the ideal MHD regime prevails, this flux can be 
written as 



' MHD 



= vB 2 - B(B • v) ~ vB 2 + B\v 8 B e 



(22) 



which shows that the energy will escape from the 
disk along flow streamlines and magnetic field lines 
(which are both directed upwards in the corona). 

For our MHD simulations presented in the previous sec- 
tion, we can define and follow the evolution of the follow- 
ing different powers. The energy liberated by accretion 



where 
Pmec = 



Plib = Pmec + Pent + Pmhd , 



Si 



dS/ • pv [ — + $ G 



Se 



Pi 



ENT : 



dS E ■ pv \— + $ G 

H« ( 1 P 

s E V7-lp 



P 



MHD : 



dS/ ■ (E x B) 



Si 



Se 



(23) 

dS E • (E x B), 



where the surfaces Si and Se correspond to the in- 
ner and outer vertical cut through the accretion disk 
(dSj = -2'KR I dZe R with -H < Z < H and dS E = 
2irREalZe.B with —H < Z < H). In standard disk 
models, the azimuthal velocity of matter is close to 
Keplerian so that the energy liberated by accretion is 

Plib ~ GM*Ma,i/2Ri where M a ,i is the accretion rate 
evaluated at the inner radius. In ADAF models, the en- 
thalpy contribution Pent is not negligible and since it 
is negative, the liberated energy (responsible for the disk 
luminosity in that particular case) is reduced. The con- 
tribution of the radial Poynting flux is more difficult to 
estimate inside the disk but we will show that it is a small 
contribution to the liberated energy. 

In the present study, the occurrence of a magnetized 
jet outflow requires that we also evaluate the energy flux 
through the disk corona. Again we define a jet power 
Pjet such that 



Pjet = Pmec, J + Pent,j + Pmhdj 
Pmec,j = J J dS surf • pv + $ G 

7 P 



Pi 



ENT,. I z 



MHD, J z 



dS surf • pv 
s mf V7 - 1 P 



dS 



surf 



(E x B) 



(24) 



where dS sur / = 2-kRoIRgz with Ri < R < Re- In 
the disk corona, the vertical velocity of matter can be 
large due to magnetic acceleration and will produce a 
significant vertical energy flux. If the structure reaches 
a near-steady state, our energy balance would read 
Pmec + Pent + Pmhd = Pjet- If neither a magnetic 
field nor a jet is present, this conservation gives the usual 
ADAF-like relation Pmec — —Pent so that Plib — 
and the disk will thus be very under- luminous. 

4.2. Under-luminous disks producing jets 

We display in Fig. [S] the temporal evolution of all pre- 
viously defined powers that characterize the energetics of 
the structure. First of all, the upper left plot represents 
the temporal evolution of the energy liberated by the ac- 
cretion of matter Plib, this power being normalized to 
the "standard" accretion power expected for a thin disk 
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Fig. 6. — Same simulation than the last snapshot of Fig^but with a square simulation domain and "force-free" boundary for the toroidal 
magnetic field (J p X B p = 0). The collimation of the jet is still occurring and the jet radius is the same than in the previous calculation. The 
dotted-dashed line is the contour of the fast-magnetosonic surface where the poloi'dal velocity of the flow equals the local fast-magnetosonic 
speed. 
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Fig. 7. — Plots of radial forces acting on the jet at Z=60 and at T=30. The total force is very small and ensures a good collimation 
for the jet. The radial equilibrium of the jet is achieved thanks to a thermal pressure gradient and centrifugal force balancing the radial 
magnetic force. This equilibrium where thermal pressure is of the same order than magnetic pressure is the signature of a "hot" jet. 



P a = GM*Mao/2Ri where Mao is the inner accretion 
rate at T = 0. This power rapidly increases in early 
stages of the evolution and then reaches a plateau where 
it remains almost constant at about 10 time its standard 
value (mainly controled by the inner accretion rate, see 
discussion in section ETT|) . In the three following plots, 



we see the relative amplitude of each contribution to this 
power. The mechanical power Pmec is the dominant 
contribution while advected enthalpy Pent and Poynt- 
ing flux Pmhd are small and of the same order, with 
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Fig. 8. — Temporal evolution of powers defining the complete energy budget as presented in Section|I] The accretion liberated power is 
non-vanishing in our structure unlike ADAF-type flow. Indeed, disk enthalpy and magnetic energy fluxes are small compared to mechanical 
power which is mainly sent into the jet. This structure is completely consistent with an under-luminous disk launching bright jets. 



ratios as 



-p, 



ENT 



0.27 



-Pi 



MHD 



0.1 



(25) 



Pmec Pmec 
Note that these two contributions are negative, meaning 
that they both advect a portion of the energy through 
the inner radius and reduce the available disk power. 

Nevertheless, these two last contributions are small 
compared to the mechanical energy released by accre- 
tion so that Plib — 0-63Pmec- This configuration is 
completely different from a usual ADAF, since there one 
would expect to have Plib "C Pmec- The jet flow is 
completely responsible for this seeming discrepancy. In- 
deed, as shown in the bottom-left plot of Fig. [SJ the 
power going into the jet (evaluated at Z = 25) is of 
the same order than Plib which leads to an energy 
budget which is indeed close to a stationary state since 
Plib = Pmec + Pent + Pmhd ^ Pjet- The jet power 
components Pmec, J, Pent.j and Pmhd.j are shown in 
the bottom-right panel: the mechanical energy flux is 
again the largest one but other components are not neg- 
ligible, in particular the jet MHD Poynting flux. This 
means that in the corona, the magnetic field is still carry- 
ing energy and consequently that acceleration of jet mat- 



ter will still occur beyond this surface located at Z = 25. 
Note also that the enthalpy flux at the base of the jet 
is a significant fraction of the jet power which confirms 
our previous statements about the "hot" nature of the 
jet. As a final illustration of the energy transfer that 
is dominated by the jet, we show in Fig. [5]both poloi'dal 
streamlines and poloi'dal MHD Poynting flux streamlines. 
We clearly see that both accretion and ejection occurs 
and that the main part of the Poynting flux is feeding 
the jet structure. Moreover, the disk remains thin (as 
seen in the flow streamlines), although Joule heating oc- 
curs within the disk. Here again, the role of the jet is 
primordial since it enables the accretion disk to transfer 
the accretion energy into jet power rather than in disk 
enthalpy (temperature) that would thicken the disk. 

5. SUMMARY AND OUTLOOK 

In this paper, we presented time-dependent MHD 
simulations of an axisymmetric, resistive accretion disk 
threaded by a bipolar magnetic field. The computation 
evolves the full set of MHD equations, including an en- 
ergy equation with spatio-temporally varying resistivity 
triggered inside the disk only (set with an a-type pre- 
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Fig. 9. — Left: poloi'dal flow streamlines and temperature. Both accretion and ejection regimes occur. The heated plasma (Ohmic 
heating and flow dynamics) in the disk give birth to a "hot" jet where temperature isocontours coincide with flow streamlines. The disk 
aspect ratio is small thanks to a predominant energy transport achieved by jet launching rather than by enthalpy creation. Right: poloi'dal 
streamlines of the MHD Poynting flux. The main part of the MHD energy flux is sent into the jet and sustains powerful hot jets. 



scription). We assumed that both the viscous torque 
and radiative losses were negligible, so that we can deal 
with radiatively inefficient, magnetized accretion disks. 
Initial conditions were carefully designed to fulfill the 
requirement of jet production, i.e. starting from a near- 
cquipartition accretion disk. Our computations bring im- 
portant new results: 

• The present simulations are the first simulations to 
achieve a near stationary state consisting of a mag- 
netized accretion disk launching trans- Alfvenic up 
to super-fastmagnetosonic, collimated, "hot" jets 
where thermal pressure balances magnetic effects. 
The initial configuration is general and valid for all 
sub-Eddington non-relativistic disks (modifications 
might occur in the very close vicinity of a black hole 
R < lORs) since it is assuming a Keplerian accre- 
tion disk in a vertical hydrostatic equilibrium and 
threaded by a bipolar magnetic field. This result 
differs from our previous study (CK02) by the in- 
clusion of a realistic energy equation replacing the 
simplifying polytropic relation. The obtained jet 
displays a temperature higher than the disk one, 
and is consistent with the creation of a hot corona. 
In this jet, we have shown that a significant part 
(~ 10%) of the jet power is enthalpy driven, and 
the thermal pressure is of the order of the mag- 



netic pressure. Moreover, the present simulations 
started from an initial magnetic configuration with 
bent magnetic field lines and therefore prove unam- 
biguously that the obtained collimation of the jet 
is self-consistently built up. We have, in addition, 
completely quantified the action of the magneto- 
centrifugal acceleration of matter, by demonstrat- 
ing the force balance along field lines as well as 
radially, showing the butterfly current circuit and 
fully diagnosing the energy transfer. 

• Evaluating energy fluxes through inner and outer 
radius as well as through the disk surface, we pre- 
sented the entire energy budget of the accretion- 
ejection structure. Unlike ADAF-like flows, we 
show that the main energy transport is achieved 
by the jet which channels the major part of the en- 
ergy released by accretion. This kind of flow can 
lead to under-luminous accretion disks supporting 
bright jets (where the assumption of non-radiative 
plasma can no longer b e valid), as for instance M87 
l|DiMatteo et al]l2f)03K In our simulations, the lu- 
minosity of the jet can be up to the energy liberated 
by the accretion Pjet ~ GM*Ma/2Ri, where M„ 
is the mass of the central object and Ri the inner 
radius of the jet. Therefore, these simulations are 
a complete model of the generic "Advection Dom- 
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inated In flow-Outflow Solutions" (AD IOS) postu- 
lated bv lBlandford fc Begelmanl l)1999fl . where the 
authors suggested that outflows might carry a part 
of the total disk energy. The accretion disk with its 
associated hot jet is no longer subject to strong en- 
thalpy creation as in ADAF flow. This also implies 
that this kind of accretion disk can remain thin, 
even if some local heating (like the Ohmic heating 
included in the simulation) is occuring inside. Fi- 
nally, these simulations adequately apply to both 
YSO and AGN type systems. 

Forthcoming work will need to deal with alleviat- 
ing the few simplifying assumptions in this work, the 
most prominent being the assumption of axisymme- 
try. This precludes our simulations from effects due 
to non-ax isymmetric instabilitie s that might perturb 
the flow llKi m fc Ostrikerl I2 000D or the equipartition 
disk ijKeppens. Casse fc Goedbloedll200^ . We will be 
able to check the stability of our simulations by stepping 
into a three-dimensional framework. Another simplifying 
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assumption was to suppose that the viscous torque is very 
small compared to the magnetic torque (as well as viscous 
heating small compared to Joule heating). The inclusion 
of viscosity in our simulations can lead to different kinds 
of flow where enthalpy creation can be enhanced and the 
jet luminosity correspondingly decreased. The transfer of 
angular momemtum will then be more complicated (both 
radial by viscosity and vertical by the jet), but will un- 
doubtedly enrich momentum, angular momentum, and 
energy transport phenomena within the accretion disk- 
jet system. 
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